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We discuss apparent paradoxes regarding the location of the zeros of the partition function in the 
complex /3 plane (Fisher's zeros) of a pure SU{2) lattice gauge theory in 4 dimensions. We propose 
a new criterion to draw the region of the complex /3 plane where reweighting methods can be trusted 
when the density of states is almost but not exactly Gaussian. We propose new methods to infer 
the existence of zeros outside of this region. We demonstrate the reliability of these proposals with 
quasi Gaussian Monte Carlo distributions where the locations of the zeros can be calculated by 
independent numerical methods. The results are presented in such way that the methods can be 
applied for general lattice models. Applications to specific lattice models will be discussed in a 
separate publication. 

PACS numbers: 11.15.-q, 11. 15. Ha, 11.15.Mo, 12.38.Cy 



I. INTRODUCTION 

The locations of the zeros of the partition function in 
the complex /3 plane contain important information re- 
garding the behavior of lattice models and statistical sys- 
tems at small and large (3 but also about possible phase 
transition or crossover behavior at intermediate values of 
f3. In the following, f3 can either be 2N / for a lattice 
gauge theory or —J/kT for an Ising model or its sigma 
model generalizations. In the following, we will use the 
notation S for the sum over all the plaquettes. For the 
Ising case, this quantity should be replaced by the sum 
of the product of neighbor spins over all the links. 

The importance of the zeros of the partition in the (3 
plane, not to be confused with the Lee- Yang zeros in the 
complex magnetic field plane, was first pointed out by 
M. Fisher for the two-dimensional Ising model and are 
called "Fisher's zeros" [J. The relationship between the 
location of the zeros and the order of the phase transition 
has been discussed in Refs. 0, 

In practice, one can calculate the partition function in 
the complex /3 plane up to an overall constant by using 
the reweighting method with respect to a real value /3o 
that is not a zero of the partition function. The partition 
function at /3 = /3o + A/3 is proportional to the average 
of exp(— A/?^) calculated at (3o: 



Zipo + A/3) = Z{po) < exp(-A/3^) >p. 



(1) 



As S has a distribution peaked near a value of the order 
of the number of sites, Z{f3o + A/3) has rapid oscillations 
when the imaginary part of Af3 is increased and it is 
convenient to subtract < S > from S. The main quantity 
that we will calculate is 



< exp [-Al3{S- < S >f)„)] >i3o 

= exp [A/3 < S >p,]Z{l3o + A/3)/Z(/3o) 



(2) 



It has the same complex zeros as Z{Po + A/3). 

In lattice gauge theory (LGT), the Fisher's zeros of 
pure (no quarks) SU{3) gauge theories with a Wilson 
action on 2 x and 4 x lattices were studied using 
Monte Carlo (MC) configurations to perform the average, 
in order to locate the finite temperature phase transition 
[3, H, Q. At zero temperature, no phase transition is 
expected between small and large /? and consequently, 
the complex zeros should stay away from the real axis in 
the infinite volume limit. 

Recently, it has been suggested 7J that a pair of com- 
plex conjugated zeros very close to the real axis could 
be responsible for the temporary behavior of the pertur- 
bative series for the average plaquette indicating a finite 
radius of convergence (power growth) rather than the ex- 
pected factorial growth. 

The existence of a zero near /3 — 5.55 + i0.12 on a 
4^^ lattice d, 0| was confirmed in Ref. @. However, it 
seemed that for a 8"* lattice, the zeros would be far away 
from the region of confidence defined in Ref. @. Using 
the same criterion, it is also not clear that the estimation 
[i| of a SU{2) zero near /3 = 2.23 + i0.155 on a 4^ lattice 
is reliable. These issues are reviewed in Sec. [Ill 

An important feature of the SU (2) distributions of ac- 
tions is that they are almost Gaussian. The distribution 
can be expressed as the product of a density of states 
n{S) (also called spectral density) and the Boltzmann 
weight exp(— /3S'). Consequently, a quasi-Gaussian distri- 
bution corresponds to a quasi-Gaussian density of states. 
In the Gaussian limit, the partition function has no zeros 
6]. 

The main goal of the article is to provide new methods 
to find the zeros of the partition function in terms of the 
small deviations from the Gaussian behavior. In order to 
be able to apply the methods to generic lattice models, 
we first introduce natural units (Sec. IIIip where the mean 
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FIG. 1: Distribution of the 1,600,000 values of S in an his- 
togram with 100 bins for an SU{2) pure gauge theory on a 4'' 
lattice at /? = 2.225 . 



and the variance are reduced to and 1. It is important 
to understand that for lattice models, the variance of S 
increases with the volume like V^^^ and that when we 
work in natural units, we need to go to larger imaginary 
part when the volume increases in order to study the 
behavior near a fixed value of f3. To make this concrete, 
in 4 dimensions, if we double the linear size of the lattice, 
the radius of confidence shrinks by a factor 4. 

The rest of the paper is organized as follows. The 
Gaussian case and the notion of circle of confidence are 
reviewed in Sec. |TVl Possible causes for underestima- 
tion of the fluctuations are discussed in Sec. |Vl A new 
criterion designed to complement the notion of circle of 
confidence for non-Gaussian distributions is proposed in 
Sec. |VT] and applied to two examples in Sec. |Vll The 
conjecture that the real part of the zero is located at a 
value of P where the variance of S is maximal is verified 
for these examples in Sec. I Villi Fitting methods to lo- 
cate the zeros in region not accessible by MC methods are 
proposed in Sec. llXl and proved successful and consistent 
for the two examples considered. 



II. THE CASE OF SU{2) 

For SU{2) LGT in 4 dimensions, action distributions 
are nearly Gaussian. This is illustrated in Fig. [TJ How- 
ever small departures from the Gaussian behavior can be 
observed. Groups of neighbor bins move above and below 
the Gaussian curve at a level that cannot be explained by 
statistical fluctuations. It is well-known Q that Gaussian 
distributions have no complex zeros. Consequently, the 
small departures should be the cause of possible zeros. 
In Ref. 9|, a zero was reported near P — 2.23 -I- i0.155 
using a reweighting of a MC distribution obtained for 
Pq = 2.225 and 39,000 configurations. This result was 
obtained by locating the intersection of the zeros of the 
real and imaginary part. An error of 0.01 on the real and 



imaginary part was quoted. 

Later, a criterion was proposed @ to determine a re- 
gion of confidence for the zeros. The criterion is based 
on the Gaussian approximation and can be summarized 
as follows. We first define 



^< > ~ < S 



(3) 



Using Gaussian integration it can be shown that for 
Nconf. independent configurations, the fluctuation in 
exp(— A/3(S'— < S >)) become of the same size as the 
average for 



|A/3|2 < ln(7Veo„/.)/4 



(4) 



This defines a circle of confidence with radius 
^ln(iVcon/.)/o'|. For the calculation of Ref. [9] , Nconf = 
39, 000 and Og = 580, we obtain a radius of confidence of 
0.135. Consequently an imaginary part of 0.155 is defi- 
nitely outside the circle of confidence. This could either 
mean that there are larger errors than the ones quoted 
or that the Gaussian criterion does not apply precisely. 
We have attempted to reproduce the result of Ref. 

using the heatbath method of Kennedy and Pendle- 
ton [l^ to generate 1,600,000 configurations. We have 
then used various methods to produce 39,000 decorre- 
lated configurations. As explained in the Introduction, in 
order to eliminate rapid oscillations we subtracted < > 
in the exponent. This does not affect the zeros of the 
partition function. In Ref. this subtraction was not 
performed and the zero level curves of the real part or 
the imaginary part appear with an approximate spacing 

01 tt/ < S >~ 0.005 in the imaginary direction. 

_ We have not found clear evidence for the result of Ref. 
9] . In Fig. [5] we show the zeros of the real and imagi- 
nary part obtained from two independent bootstrap pro- 
cedures where the configurations are picked at random 
with possible repetitions from the large original set of 
1,600,000. 

In both graphs, no zeros are found within the error 
bars of the Ref. y\. In addition, zeros are found slightly 
within the circle of confidence and at different locations. 
For comparison, the figures show the location of the zeros 
of the real and the imaginary part that would be seen if 
the distribution was a Gaussian with the same mean and 
variance than the sample. These zeros appear on non- 
intersecting hyperbolas ^ given by the formula 



(5) 



with odd (even) n standing for the zero of the real (imag- 
inary) part. One robust feature found in the collection 
of figures produced with independent bootstraps is the 
way the MC data departs from these hyperbolas. On the 
far right the MC data is above the hyperbolas and then 
goes below as we move toward the center. On the left, 
the MC data follows more closely the Gaussian hyperbo- 
las. In summary, the MC data seems slightly tilted to 
the left. Another robust feature which seems correlated 
with this left tilt is that the complex zeros within the 
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FIG. 2: Zeros of the real (circles) and imaginary (crosses) 
part for 39,000 configurations of a SU{2) gauge theory on a 
4'* lattice with /3o ~ 2.225. The solid lines are the circle of 
confidence and the hyperbolas of the Gaussian approximation 
described in the text. 



models, for instance scalar models or spin models where 
the typical scales are completely different. We define a 
reduced form: 



Sred. = (S- < S >)/as 



(6) 



As < 5* > is already subtracted in the average consid- 
ered in Eq. ([2]), we only need to rescale A/3 in order to 
eliminate the variance: 



f3r 



A(3as . 



(7) 



For instance, in Fig. [21 the scale of the figure is about 
0.16. In the reduced coordinates, the size of the figure 
would be about 3.85 and could be compared with corre- 
sponding graphs for other models. 

In order to make some equations more readable we 
introduce the following notations: 



Pred. — 


X + iy 


/ ^ 


< ex.p{-l3red.Sred.) > 


R EE 


Ref 
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< iR.eC^p{~Pred.Sred.) 
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"Im — 


< {ImCXp{-Pred.Sred.) 











Gaussian circle of confidence appear on the left side of 
the figure. 

At this point, it is not clear to us that with the ex- 
isting information it is possible to draw solid conclusions 
regarding the existence, the location or the absence of 
complex zeros in the portion of the complex /3 plane con- 
sidered. Consequently, we decided to study examples of 
quasi Gaussian distributions where the location of the ze- 
ros can be determined by independent numerical meth- 
ods. Another remark, understanding the reweighting in 
the complex plane may also be important for problems 
with a chemical potential. 



III. NATURAL UNITS AND NOTATIONS 

In the following we will consider MC distributions gen- 
erated for the purpose of testing general methods that 
can be applied for a broad range of models. In order to 
translate the results presented hereafter in a form that is 
useful for specific models, we will redefine the origin and 
the scale of S. 

We have seen in the previous section, that the distribu- 
tion of Wilson action S is approximately Gaussian. Given 
this distribution, it is always possible to take every value 
of S', subtract the mean and then divide by the variance. 
We then obtain a distribution which is approximately a 
normal Gaussian. The advantage of using this normal 
form is that the results that follow can be used for other 



It is understood that f,R... etc are all functions of f3red.- 
In the following, unless specified, we will work mostly 
with Sred. and Pred. and the subscripts "red." will be 
dropped. 

IV. THE GAUSSIAN CASE REVISITED 

A. Data generation 

We used the Metropolis algorithm in order to generate 
values of S with a normal distribution: 

PiS) = (2^)-i/2 exp(-5V2) • (8) 

The upgrade S S + SS were performed with SS uni- 
formly distributed in an interval [—R,R]. The probabil- 
ity of accepting the change (acceptance) is 

dS / d(5S'exp(-S'V2) 

-oo J-R 

[0{dS)9{2S + 6S) + 9{-SS)9{-2S - 6S)] . (9) 

We have generated 1,600,000 values of 5 for i? = 1. The 
average of this set of values was 0.00087 and the accep- 
tance was 0.8047 in good agreement with Eq. ([9]) which 
predicts 0.804583. The values were sorted in 100 bins of 
equal width as shown in Fig. [3] The departures from 
the Gaussian behavior are difficult to see. Before analyz- 
ing these departures, we will first address the question of 
autocorrelations of successive configurations. 
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FIG. 3: Distribution of the 1,600,000 values of S in 100 bins 
(above) and the corresponding (below). 



FIG. 4: ln(|C; I) versus / for the original set of 1,600,000 values. 



B. Autocorrelations 

The autocorrelations are defined as usual by the for- 
mula 

n 

Cl = (1/n) J2 (Sm- < S >)iSrn+l- < S >) , (10) 
m— 1 

where the subscript refers to the configuration number. 
The numerical values are shown on a logarithmic scale 
in Fig. m The linear part of the graph stabilizes when 
n is large enough (here about 200,000). The slope of 
the linear part is about -0.12 which means that the cor- 
relation time is about 8. In order to select reasonably 
decorrelated values it is necessary to select values that 
are several correlation times apart. In the following, we 
constructed subsamples of size 40,000 (about the same 
size as in [S]), obtained by selecting one in every 40 val- 
ues from the original sequence. As shown in Fig. [5l 
the subset show no apparent autocorrelations. Note that 
results similar to the ones obtained in the rest of this 
section can be obtained by using the bootstrap method 
instead of the skipping method. 



C. Histogram analysis for decorrelated data 

The analog of Fig. [3]for the subsample described above 
is shown in Fig. [S] More information can be gathered by 
displaying the differences in units of the expected fluctu- 
ations ri for each bin, namely 



Gaussian MC 



= {N, - NPi)/./NP, 



(11) 



where N is the total number of values (40,000), P^, the 
expected probability to be in the i-th bin, the probability 
density (for a normal Gaussian) in the middle of the bin 
multiplied by the bin width, and Ni the number of values 
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FIG. 5: In(|Ci]) versus I for the subset 5*1, 541, . . . (above). 



in the i-th bin. If the N processes were independent, 
we would have a binomial distribution and the average 
number in the i-th bin would be NPi with a variance 
af = NP,{1 - Pi) ~ NPi, since the Pi are smafl. We 
expect I Til ~ 1. For the subset obtained with the skipping 
method J^rf ~ 62 (for 100 bins). 

It is important to notice that the show almost no 
coherent structure. The difference tend to have opposite 
signs and no simple pattern seems to be visible except 
for the tails of the distributions. Coherent negative dif- 
ference are seen for \S\ > 4. This feature has a very sim- 
ple explanation: there are almost no values for \S\ > 4 
and we see merely minus the square root of the Gaussian 
distribution. Given that we only only use 40,000 configu- 
rations, large values of \S\ corresponds to smaU statistics. 
The probability of having \S\ > 4 is about 6 x 10~^ for 
a normal distribution. Consequently, we would expect 2 
values in this range in a typical MC run described above. 
In Fig. [51 one value near 4.3 produces a positive spike. 
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FIG. 6: Distribution of the 40,000 values of Si, 5*41, 
100 bins (above) and the corresponding (below). 




FIG. 7: Zeros of the real (circles) and imaginary (crosses) 
part for 40,000 Gaussian configurations. The solid lines are 
the circle of confidence and the hyperbolas of the Gaussian 
approximation described in the text. 



V. ERROR ANALYSIS FOR THE GAUSSIAN 
MODEL 

As already shown in Ref. Q, all the quantities of in- 
terest can be calculated using standard methods of Gaus- 
sian integration. Remembering the notations introduced 
in Eq. The average and variance of the real and 

imaginary part can be expressed as: 



/ 
R 



J2 



cxp((l/2)(a;2-2;2))sin(xy) 
(1/2) exp(2a;^)(exp(-22/2) cos(4xy) + 1) 
~exp(a;^ - y^)cos^(a;?/) (12) 
(1/2) exp(2a;2)(- exp{-2y^) cos(4a;y) + 1) 
— exp(x^ — y^) sin^(a;y) 



D. Zeros of < exp(-/3S) > 



Having made a reasonable effort to produce samples 
of 40,000 values of 5* that are distributed like a normal 
Gaussian, we now consider the zeros. The zeros of the 
real and imaginary part of < exp(— /JS") > are displayed 
in Fig. [7] for 40,000 configurations. For comparison, we 
also show the hyperbolas Im(3Re(3 = mr that should be 
followed by these level curves and the circle of confidence 
{Re(3f + [Imfif = ln(40,000) ~ (3.26)2. We see that 
as we follow the level curves from the inside of the circle 
toward the boundary, the agreement slowly deteriorates. 
The apparent complex zero closest to the origin is near 
(3 ~ 1.5 -I- z3.0 and is slightly outside the circle of confi- 
dence. Fig. [7] is thus consistent with the known absence 
of zeros in the Gaussian case. Note that it is not impossi- 
ble to find random samples of 40,000 configurations with 
complex zeros slightly inside the circle of confidence, but 
this is rather uncommon. By design [(|, this should occur 
in 16 percent of the cases in average. 



We have compared these analytical formulas with MC 
results for x — Re(3— 1, 2 and 3. The errors bars on the 
MC data point have been estimated by using the MC 
estimator of the variance divided by ^40, 000. 

The results for the real part at fixed a; =1 are displayed 
in Fig. m One sees that the there is a very good agree- 
ment until the function becomes smaller than the error 
bars. This occurs near y = 3 which is approximately 
where the crossing of the circle of confidence occurs. A 
more detailed investigation also shows that for this value 
of X, the MC estimation of the variance is in good agree- 
ment with Eq. (jl2p within the circle of confidence. 

Larger discrepancies and larger error bars can be seen 
in Fig. [9] for x = 2. At some places the analytical ex- 
pression is almost outside the MC error bars. Also, the 
MC estimate of the variance is now significantly lower 
than the analytical one as shown in Fig. [TOl In other 
words, the errors bars of Fig. [5] should be larger. All 
the features observed at a; = 2 are dramatically ampli- 
fied for a; = 3. The analytical expressions is now mostly 
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FIG. 8; R — Re < exp(— /SS* > versus y — Imp at fixed 
X — Rep— 1. The continuous line is Eq. (|12|l . the circles are 
MC results. The lower figure is a magnified part of the upper 
figure. 
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FIG. 9: R — Re < exp(— /JS* > versus y — Imf3 at fixed 
X = Rel3= 2. The continuous line is Eq. p2|l . the circles are 
MC results. 



outside of the MC error bars in Fig. [TTJ A more detailed 
study shows that the error bars are severely underesti- 
mated (sometimes by an order of magnitude). The dis- 
crepancies observed in Fig. [Tl]have a simple explanation. 
When Rej3 becomes large, the configurations with large 
negative values of S become very important. We have 
seen in Sec. IIV CI that no configurations with S < —4 
were present (while we should typically expect one). Fig. 
[T^ shows that the contribution of the tail S < —4 fits 




FIG. 10: cr|j versus y = Imf3 at fixed x = Re/3^ 2. The 
continuous line is Eq. (fT2)> . the circles are MC results. 
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= Re < exp{—/3S > versus y — 
The continuous line is Eq. (fT2)). 
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very precisely the discrepancy in the average value. The 
absence of configurations with S < —4 also explains the 
underestimation of the error bars. 

In summary, when Re(3 is not too large, the MC results 
are in good agreement with Eq. within the circle of 
confidence. When \Re/3\ becomes too large, we probe the 
tails of the S distribution where we have low statistics 
and obtain unreliable results even within the circle of 
confidence. 

A solution to the problem of low statistics consists in 
introducing a bias exp(— 65) in order to probe better the 
tails of the distribution. In terms of the original prob- 
lems discussed in the introduction this consists in using 
reweighting for different but nearby values of (3o. This 
amounts of "patching" the various data as explained in 
Refs. [l,|Tl[. The question will be discussed in Sec. I Villi 
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FIG. 12: Discrepancy between Eq ^ and MC in Fig. [TT] 
(solid square) and analytical contribution of 5 < —4 to i? = 
i?e < exp(-/35' >. 



VI. A NEW CRITERION FOR THE REGION OF 
CONFIDENCE 

We have seen in the previous sections on the Gaussian 
distribution, that provided that i?e/3 is not too large, the 
MC results are reliable within the circle of confidence 
and that it is uncommon to produce complex zeros in 
this region. On the other hand for SU{2), complex zeros 
within the circle of confidence are very common. Since in 
this case the distribution is not exactly Gaussian, com- 
plex zeros are possible but the changes observed in Fig. [2] 
suggest that the complex zeros seen in the circle of con- 
fidence (calculated in the Gaussian approximation) are 
numerical artifacts. In general, we may obtain false com- 
plex zeros if the zero level curves of the real and imagi- 
nary part move too far from their exact trajectories. In 
this section, we propose a criterion to decide when it is 
the case. 

The Gaussian circle of confidence in the complex /3 
plane is defined by the condition 

^/(/3)/\/^W < l/(/3)l (13) 

This works very well in the Gaussian case because /(/3) 
is never zero. On the other hand for a non-Gaussian dis- 
tributions with complex zeros, this criteria will exclude 
a neighborhood of the zero. 

We would like to know how the MC level curves move 
as a consequence of the uncertainty in /. A possible 
method to search for the zeros of say the real part con- 
sists in following this quantity in a particular direction 
and look for the change of sign. The error on the loca- 
tion of the zero is then the error on the function divided 
by the slope. This error depends on the direction. It 
can be minimized by moving in the normal direction to 
the level curve where the slope is maximal. A complex 
generalization of this estimate is the ratio tT/(/3)/|/'(/3)|. 
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FIG. 13: Boundary of the confidence region for d — 0.3 (cir- 
cles, upper set), 0.2 (crosses, middle set) and 0.1 (boxes, lower 
set), compared to the Gaussian circle of confidence, all for 
40,000 configurations. 

This defines direction independent circles of uncertainty 
around every point on the complex plane. We propose 
to consider the alternative regions of confidence defined 
by the condition. 



In order to be useful d should be a fraction of the typ- 
ical distance between zero level curves of the real and 
imaginary part. Note that for the Gaussian distribution, 
/g(/5) = exp[(l/2)/32] and consequently = /3/g. If we 
impose that the the two conditions should be equivalent 
on the original circle of confidence Eq. (|13p . we obtain 
d ~ 0.3 for N = 40, 000. We have checked (Fig. [TS]) that 
this condition produces a region of confidence very close 
to the Gaussian circle of confidence for \/3\ < 2. We have 
seen in Sec. |Vl that for \f3\ > 2 we usually underestimate 
the variance and the region of confidence becomes larger 
than it should be. 

Note that Eq. ([TH)) guarantees at 84 percent confi- 
dence level that we do not have artificial zeros. For the 
purpose of establishing lower bounds, one might want 
have a higher confidence level. The same remarks apply 
for the new criterion. Empirically, we found that d — 0.2 
removes an upper sliver of the Gaussian confidence circle 
where most of the occasional artificial zeros appear. This 
is illustrated in Fig. [131 

A practical remark. In order to calculate |/'(/3)|, we 
can make use of the fact that / is an analytical function. 
This is obviously true in the Gaussian case. This is also 
obviously true for any average over a finite number of 
configurations. Consequently, we may use the Cauchy- 
Riemann conditions. 

Assuming analyticity, one can calculate the derivative 
in an arbitrary complex direction. For instance, if we 
search for zeros at fixed values of i?e/3, it is convenient 
to use the derivative in the imaginary direction. In this 
particular case, we have 

l/'l = y/idR/dyf + {dl/dyf (15) 



8 



VII. EXAMPLES OF QUASI-GAUSSIAN 
DISTRIBUTIONS 

We now consider normal Gaussian distribution with 
small perturbations of order 3 and 4: 

P{S) cx exp(-(l/2)S'2 - XsS^ - U"^) (16) 

Note that due to the non-Gaussian corrections, < S >^ 
and the variance is usually not 1. For this reason, we use 
the tilted symbol to emphasize that the transformation 
Eq. ([6]) needs to be used in order to have the natural 
units advocated in Sec. IIIII After this transformation is 
performed we get 

P(S) cx P{<S> +(JsS) 

cx exp(-Ai5 - AaS'^ - AaS"^ - A4S"') (17) 

where the untilded symbol corresponds to the reduced 
units (but the "red." subscript is dropped). 

A. Example 1: A3 = 0.1, A4 = 0.01 

We first consider the case A3 = 0.1, A4 = 0.01. It has 
been chosen in such a way that we have zeros inside and 
outside the Gaussian region of confidence. It is possible 
to find the zero level curves of the real and imaginary 
parts by using numerical integration. We refer to these 
curves and to the complex zeros where they intersect as 
"accurate" . The accurate zero level curves are compared 
with the MC ones in Fig. [TH The accurate complex zeros 
are located at f3 = 1.3735 + il.7104 and /3 = 1.3735 + 
«2.9478. As far as we can tell the two real parts are the 
same. For further reference, the variance of the original 
variable S is 1.0988 and the location of the real part 
of the complex zeros is 1.2500 in unrescaled units. The 
MC zeros closest to these accurate zeros are located near 
1.3+11.6 (close to the first one) and 0.9+12.8 (not so close 
to the second one). We have considered several values of 
d to define a region of confidence. It clear that it should 
contain the lowest zero which is a good approximation. 
On the other hand, the second zero is produced by curves 
that depart significantly from the accurate ones. A value 
oi d = 0.12 defines a confidence region that contains the 
lowest MC zero but unambiguously excludes the other 
zero (see Fig. [Ti)) . Note that there is a small island 
of exclusion not so far from the lowest zero. This is in 
agreement with larger departures where the curves turn 
over. When d is increased, the upper boundary moves 
up and the island shrinks. For d > 0.2, the confidence 
region includes the second zero. 

B. Example 2: A3 = 0.01, A4 = 0.002 

Our second example is A3 = 0.01, A4 — 0.002. The 
perturbation is much smaller and the first accurate zero 



Example 1 ; (a/\f\) <0,12 



o Re. 
X Im. 




FIG. 14: Zeros of the real (circles) and imaginary (crosses) 
part for 40,000 configurations corresponding to the first ex- 
ample. The small dots are the accurate values for the real 
(gray, green on-line) and imaginary (black, blue online) parts. 
The exclusion region boundary for d = 0.12 is represented by 
boxes (red online). The solid line is the circle of confidence of 
the Gaussian approximation. 



Example 2 ; (tj/|f1) <0.15 




Se. |S 

FIG. 15: Zeros of the real (circles) and imaginary (crosses) 
part for 40,000 configurations corresponding to the second 
example. The small dots are the accurate values for the real 
(gray, green on-line) and imaginary (black, blue online) parts. 
The exclusion region boundary for d = 0.15 is represented by 
boxes (red online). The solid line is the circle of confidence of 
the Gaussian approximation. 

is at /3 = 1.207 + i5.241, far outside the Gaussian circle 
of confidence. All the zeros seen in Fig. [15] are artifacts. 
The value d = 0.15 excludes all them with a good margin 
of safety. 

VIII. CONJECTURE CONCERNING THE 
REAL PART OF THE COMPLEX ZEROS 

In Sec. |Vl we mentioned the possibility of probing the 
tails of the distribution by introducing a bias. In this 
section, we will explore this possibility by introducing a 
real parameter b to favor one side of the distribution: 

Pt{S) = P{S)exp(-bS) . (18) 

In the language of the original problem it amounts to 
shift the value of /3o- In this section, we do not perform 
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Example 1 with a bias b 




FIG. 16: Second moment as a function of b. 



the change of variables as in Eq. ([5]), because we want 
to discuss the b dependence of the variance (which would 
become 1 after the rescaling) . For this reason we keep the 
tilde in the equations. If we denote the average with the 
new distribution as < ... >{,, we have the simple identity 



< 



exp(-/35) >o= C{b) < cxp(-(/3 - b)S) >b . (19) 



C{b) is a normalization factor that will be assumed non- 
zero. This assumption can be justified as long as there 
are no zeros of the partition function on the real axis (we 
remind that b is real). It is then clear that the effect of 
the bias is to shift a complex zero by -b. 

In the Gaussian case, the zeros have a mirror symmetry 
with respect to the imaginary axis. When odd perturba- 
tions are introduced, the symmetry is broken. In Figs. [HI 
and 1 151 the complex zeros are on the right and the curves 
are in some sense tilted to the right. As we increase b by 
positive values, the complex zeros move to the left. The 
third moment of the distribution, < {S— < S >& can 
be seen as a measure of the left-right asymmetry. We con- 
jecture that when the real part becomes zero, the third 
moment vanishes, which corresponds to an extremum of 
the second moment. This conjecture has been checked 
very precisely with the two examples of Sec. IVIII It is 
illustrated in Fig. [1^] for example 1. We see that the 
second moment peak near 1.25 which coincides with the 
location of the real part (before rescaling) of the complex 
zeros. 

If this conjecture remains true even approximately, it 
is telling us to look for complex zeros for real values on a 
vertical strip that includes the region(s) of the real axis 
where the second moment is maximum. 



IX. INDIRECT METHODS TO LOCATE THE 
ZEROS 

So far, we have not discussed methods to find com- 
plex zeros with a large imaginary part. This may occur 
as in Example 2 when the corrections to the Gaussian 
form are small. Under these circumstances, the third 
and fourth moments, that would vanish in the Gaussian 





MC 


Accurate 


ma 
m4 


-0.575 
0.498 


-0.559 
0.448 



TABLE I: Third and fourth moments from MC and accurate 
numerical integration in example 1. 





MC-l-Newton 


Accurate 


Ai 


-0.320406 


-0.315192 


Aa 


0.50140 


0.496627 


A3 


0.117795 


0.115924 


A4 


0.01433 


0.014577 



TABLE II: Ai from Newton's method compared to their ac- 
curate values in example 1. 



case, are possibly very small and it might be possible 
to infer the small non-Gaussian corrections from these 
small moments. For completeness, the third and fourth 
moments are defined as 



ma 
mi 



< ^red. 



> 



> 3 < S'^ed. >^ 



(20) 



We have used Newton's method to determine the un- 
known parameters Ai, . . . A4 in Eq. (|17p from the first 
four moments calculated with MC. This is the only in- 
put and it can be calculated easily for more complicated 
models. The numerical estimates of the moments are 
based on a single set of 40,000 configurations obtained 
as described above. We worked with the reduced vari- 
ables so the first two moments are fixed to and 1. The 
integrals and the derivatives in Newton's method were 
calculated numerically. In example 2, the corrections are 
small and the initial values (0,1/2,0,0) corresponding to 
the normal Gaussian is a good approximation. Newton's 
method converges rapidly for these initial values. For ex- 
ample 1, the corrections are larger and one needs to start 
closer to the solution. The main source of uncertainty 
come from errors in the determination of the moments 
from MC data. The numerical results are shown in Ta- 
bles 1 to 4. 

Using the approximate values of Ai obtained from the 
MC moments and Newton's method, we were able to 
construct approximate zero level curve for the real and 
imaginary parts the results are shown in Fig. [T71 In ex- 
ample 1, the agreement between the two sets of curves 
is very good and allows an accurate determination of the 
first two zeros. In example 2, the agreement is also very 
good except near the complex zero where the discrepancy 
is visible. The value of the complex zero obtained with 
the new method is 0.98-1-15.14 which differs from the ac- 
curate value 1.207-1-15.241. We have repeated the same 
calculation with other sets of 40,000 configurations ob- 
tained with the bootstrap method and obtained similar 
size for the discrepancies. 
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MC I 


MC II 


Accurate 


ma 

7714 


-0.0469 
-0.0408 


-0.0576 
-0.0371 


-0.0535 
-0.0353 



TABLE III: Third and fourth moments from MC I (40,000 
data points), MC II (1,600,000 data points) and accurate nu- 
merical integration in example 2. 





MC I-HNewt. 


MC II-fNewt. 


Min. 


Accurate 


Ai 


-0.0247 


-0.0303 


-0.0299 


-0.0280493 


A2 


0.488 


0.489 


0.490 


0.489364 


As 


0.00837 


0.0103 


0.0101 


0.00948691 


A4 


0.00209 


0.00206 


0.00189 


0.00192227 



TABLE IV: \i from Newton's method with the two sets of 
MC data described in the text, and from th minimization of 
compared to their accurate values in example 2. 



We have tried to improve the accuracy of the estimates. 
First, we have used the entire data (1,600,000 values) to 
determine the moments. This improved the accuracy of 
the moments (see Table III) and of the parameter Ai (Ta- 
ble IV). The complex zero was located at 1.21-1-15.08. We 
have also used a minimization of the of the histogram 
to determine the Aj. The results are very good (see Table 
IV) and produces level curves that a almost impossible 



to distinguish from the accurate ones on a graph. 

As the perturbations get smaller, the zeros get a larger 
imaginary part and the numerical integration becomes 
more difficult because of the fast oscillations of the inte- 
grand. However, it possible to use perturbative methods. 
When A3 and A4 are both zero, the problem is Gaussian, 
solvable analytically and there are no complex zeros. If 
we calculate < exp(— /35') > at first order in A3 and A4 
and divide by the Gaussian limit (which has no zeros), 
we obtain a polynomial of order 4 in /3: 



< exp(-/3S')(l - X^S^ ~ A4S"*) >G / < exp(-/3S') >g 
= QW 

= ! + ■■■- A4/3V(16A^) , 



where < • • • >g is a notation for the average with a 
Gaussian weight as in Eq. pT|) with A3 and A4 set to 
zero. After expansion, the polynomial Q{(3) has 16 terms 
which can be calculated easily by Gaussian integration. 
For the accurate values of of Table IIXI we have zeros 
of Q{P) = for /3 = 1.08 ± i4.56 which is not far from 
the correct value. 



Example 1; Zeros from fit 



O Re, 
X im. 




0.0 h , 1 ^ 1 ^ 1 ^ I , , , 1 , , , c 

-4.0 -3.0 -3,0 -1.0 0.0 1.0 2.0 3.0 4.0 
Re. (9 

Example 3; Zeros from fit 




0.0 tl ' 1 • 1 • 1 ' B ' 1 ' 1 ' 1 • C 

-4.0 -3.0 -3.0 -1.0 0.0 1.0 3.0 3.0 4.0 
Ee. /S 



FIG. 17: Zeros of the real (circles) and imaginary (crosses) 
using the approximate Ai obtained from MC moments. The 
small dots are the accurate values for the real (gray, green on- 
line) and imaginary (black, blue online) parts, for example 1 
(upper graph) and example 2 (lower graph). 



X. CONCLUSIONS 

We have build a "ladder" of methods that can be ap- 
plied for increasing values of the imaginary part in natu- 
ral units. We have refined the determination of the region 
where MC calculation can be trusted. Fitting methods 
based on cubic and quartic perturbations work for larger 
values of the imaginary part. It is clear that the paramet- 
ric form of the distribution is known in the two examples 
considered here. When the form is not known, this in- 
troduces an additional uncertainty. The zeros for the 
approximate distributions can be found using numerical 
integration provided that the imaginary part is not too 
large (about 5 to 6 in natural units). Larger values in- 
dicate small non-Gaussian terms and can be treated by 
perturbative methods. These methods are being experi- 
mented on SU{2) and SU{3) LGT. This wiU be reported 
in a subsequent publication. 
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